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It is shown that the generalization of the Navier-Stokes 
equations to a theory with N "internal state" copies of the 
velocity fields is a step in a wrong direction: the N — » oo limit 
has no physical sense and produces wrong results, whereas the 
treatment of the first order terms in 1/N is even more com- 
plicated than the initial problem of description of turbulence 
in the frame of the Navier-Stokes equation. 



Consider an incompressible velocity field u(r, t) which 
is the solution of the Navier Stokes equations 

du/dt + (u ■ V)u - vV 2 u + V P = , V-u = 0, (1) 

where v is the kinematic viscosity and p is the pressure. 
We assume that there exist appropriate boundary condi- 
tions to maintain a high Reynolds number flow. It has 
been suggested in the literature, and especially recently 
[jlj , that it is advantageous to consider a generalization of 
this equation to a situation in which there are N copies 
of the velocity field labeled by an "internal" index s: 



u(r,t) u s (r,t), s=l,2, . 



(2) 



In terms of these copies one writes the generalized equa- 
tion for incompressible fields u s : 

du s /dt + A slm (ui ■ V)u m - vV 2 u s + Vp s = . (3) 

The hope is that the statistical properties of this theory 
are simpler to elucidate than those of (Q) in the limit 
TV — > oo. We will show here that this hope is not realized 
when the proper symmetries are taken into account. 

Consider first the implication of Galilean invariance. 
Our theory has to remain invariant to the transformation 



r = r - U n t . 



Us 



h, ■ Un 



(4) 



In (^) h s is a scalar in 3-space, but a vector in the inter- 
nal state space. In case when the velocity fields u s have 
a physical sense (like in the two-fluid models of plasma) 
Galilean invariance requires h s = 1 for every s. However 
our future conclusions are independent of the particular 
choice of h s , which may be considered as a free parame- 
ter. 

Applying the transformation (Eh, to (j^) we find two 
additional terms. The time derivative leads to an ex- 
tra term which — (Uq • V)u s , whereas the nonlinear term 



leads to the extra term A s i m hi(XJo ■ V)u m ). Since these 
two terms must cancel, we find the constraint on A s i m 
which follows from the fundamental symmetry of hydro- 
dynamics - Galilean invariance: 



(•5) 



Next remind that in the known cases of simplification 
of a problem in the large N limit (see e.g. [|[D) one 
needs to have some continuous group of symmetry, the 
Lie group. Important example of the group Lie are the 
group of rotations in 3 dimensional space SO (3) and the 
group of unitary matrixes 2x2 SU(2). Elements of the 
Lie groups g may be parametrizied by a set of continu- 
ous parameters, like the Euler angles for the SO (3) group. 
Following and aiming possibility of 1/iV-simplification 
we equip the copy space with a symmetry group Lie Q. 
By analogy with the theories ||,|| the index s have to be 
considered as the label of the basis of the representations 
of Q. Denote by T(g) the operator that corresponds to 
the the element g of the group, and the transformed field 
by u s . Then 



u s = (f(g)u) s = T ss/ (g)u' s 



(6) 



where N x N matrix T ss i (g) is the representation of T(g) 
in this basis. Now we apply T(g) to the equations of 
motion (||) and demand invariance. Write for that these 
equations for u s ' and multiply T ss i (g) from the left. Then 



+ T ss ,(g)A snm (u r V)u m - jA7 2 u s + Vp 8 



0. (7) 



Next we use the fact that Q is a group and therefore the 
matrix T ss i has an inverse u s = T~}(g)u s '. Substituting 
in (0) and demanding invariance leads to the constraint 



T SS '(g)A s , l , m ,T l , l 1 (g)T m } m (g) = A slr , 



(8) 



This is the constraint on A s i m which follows from the 
fact that the set of Eqs. (^) is invariant with respect to 
the transformation (|J) of the group Q. 

Now we have two constraints (|^) and (§|) on the same 
tensor A s i m which will be considered as a restriction on 
the allowed form of the transformation (JsJ) of the group 
Q. To find this restriction we multiply the left hand side 
of (||) by hi and sum up on I. Together with it gives: 



T ss ,(g)A s , l>m ,T- l 1 (g)T-} m (g)h l = 5, 



(9) 
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This equation is simplified, by multiplying on the left by 
and on the right by T to obtain A sm i>Tu, (g)hi = 
5 sm . It appears now that if we have uncountably 
many constraints (since g is continuous) on the finite 
dimensional tensor A sm i- The only way to remove 
the over determination is to demand that (g) is g- 
independent. Because for the identical transformation 
(g = e) T ss > (e) = 5 SS > we have 

T ss ,(g) = 5 ss , . (10) 

Finally, if we apply (??) to (^) we find that the equipment 
of the copy space with a continuous symmetry group 
leaves A s i m without any additional constraint. Equation 
(^) becomes an empty identity. The way to understand 
this startling result is that the requirement of Galilean 
invariance introduced an anisotropic ray h s in the copy 
space, and there does not exist a nontrivial transforma- 
tion that leaves it invariant. In other words the require- 
ment of the (/-symmetry of the set of Eqs. (|J) itself gives 
nothing because the whole problem includes the Galilean 
invariance (Q) which contradicts to this symmetry. 

The point to understand now is that the conclusion 
that the coefficients A s i m are unconstrained eliminates 
any hoped for advantages of the 1/N expansion. The 
technical remark is that in successful applications of this 
method the one-loop diagrams have an TV weight that is 
larger than that of two or more- loop diagrams. The way 
this works in practice II is through some "coherency" 
conditions in A s i m that cause a 2-loop diagram to count 
N less than in a product of two 1-loop diagrams. How- 
ever, if A s i m is unconstrained by the choice of the sym- 
metry group, there is no loss of N factors in the 2-loop di- 
agrams compared with the 1-loop diagrams, and no sim- 
plification of the diagrammatics appears in the N — > oo 
limit. 

Yet in a recent contribution there was an attempt to 
overcome this hopeless situation by a clever trick. The 
idea of Q| is to consider the TV-dimensional space (||) 
of additional unphysical vector fields u s (r,i) in which 
Galilean invariance (Q) is not respected, coupled to one 
physical velocity field Uo(r,i) for which Galilean invari- 
ance is retained. The result of this attempt was a pre- 
diction that the scaling exponent £2 of the second order 
structure function in the N — > 00 limit attains the value of 
C2,iV->oo = 1/2 = 0.5 which differs from the experimental 
value C2,exp — 0.70 (see e.g. 0) much large than the well 
known Kolmogorov 1941 prediction C2,K4i = 2/3 ~ 0.67. 
There is an immediate reason for worry about this predic- 
tion. It has been observed many years ago by Kraichnan 
H that disregarding the Galilean symmetry and trun- 
cating the diagrammatic at a finite order produces a 1/2 
prediction for £2 ■ Is it possible that the result reported in 
jlj is essentially identical to that? We think so. In fact, 
on page 3761 of § one finds the following observations: 
"Since we have seen that the zero mode [the physical ve- 
locity Uo(r,<)] contributes negligible to the internal bonds 



of a graph when N — > 00, we may now completely ne- 
glect the zero mode in resuming these diagrams. This 
means that the theory for the N unphysical modes in 
the N — > 00 limit is decoupled from the physical velocity 
field and becomes a theory in which Galilean invariance 
has been totally discarded. Then we find u In particular, 
we shall now see that the Green's function and the double 
correlator [of the physical velocity Uo(r, t)] may be ex- 
pressed completely in terms of the Green's function and 
double correlator [of the unphysical modes u s (r,t)]. It 
follows from the possibility to neglect (for N — » 00) the 
own nonlinearity of of the physical velocity with respect 
to N nonlinear contributions of unphysical modes. Thus 
the scaling behavior of the physical velocity is totally 
determined by the unphysical modes. The conclusion is 
that the Galilean invariance disappears from the problem 
in the formulation of Ref. jp in N — > 00 limit. 

Clearly for N = we have the initial Navier- 
Stokes based formulation of the hydrodynamic turbu- 
lence. Therefore one may hope to have something reason- 
able in the next terms oi 1/N expansion when some 
tail of the Galilean invariance will recover. On the page 
3745 of § we found: 11 A real test of our approach would 
be to compute the first correction in powers of 1/N ..." . 
This corrections originate from not only two-loops dia- 
grams of unphysical fields but also from ALL ORDERS 
diagrams with respect to own nonlinearity of the physical 
velocity field. There are no 1/N simplification in the later 
series. Therefore the task to find the 1/N corrections is 
equivalent to solution the problem of the Navier-Stokes 
hydrodynamic turbulence. 

In summary, it was suggested the multicomponent gen- 
eralization of the problem of turbulence such that the 
zero-order approximation is solvable. However the next 
step occurs to be as complicated as the whole solution 
of the initial problem. It is not unexpected and hap- 
pens always when the zero-order problem has no relation 
to the initial one. In considered case this is so because 
the Galilean invariance (broken in zero-order) as stressed 
years ago by Kraichnan, is crucial for hydrodynamics. In 
short, the multicomponent extension of the velocity field 
discussed in JjJ is a step in a wrong direction which makes 
the following steps even more complicated. 
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